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The SU(n) Heisenberg model represented by exchange operators is studied by means of high- 
temperature series expansion in three dimensions, where n is an arbitrary positive integer. The 
spin-spin correlation function and its Fourier transform S(q) is derived up to 0[(/3J) 10 ] with /3 J 
being the nearest-neighbor antiferromagnetic exchange in units of temperature. The temperature 
dependence of S(q) and next-nearest-neighbor spin-spin correlation in the large n cases show that 
dominant correlation deviates from q = (71", n, 71") at low temperature, which is qualitatively similar 
to that of this model in one dimension. The Neel temperature of SU(2) case is precisely estimated 
by analyzing the divergence of S(ir,n,n). Then, we generalize n of SU(n) to a continuous variable 
and gradually increases from n = 2. We conclude that the Neel ordering disappears for n > 2. 

PACS numbers: 75.10.Jm 75.40.Cx 75.50.Ee 



I. INTRODUCTION 

It is known that properties of quantum spin systems 
tend to approach those of their corresponding classical 
spin systems as the spin magnitude increases. How- 
ever, this is not necessarily the case for a sequence of 
models in which the number of multipolar-interaction 
terms increases as the spin magnitude increases. Conse- 
quently the higher-spin systems may have stronger quan- 
tum effects in this case. In other words, such addi- 
tional terms can break a classical correspondence down 
even in high dimensions. These "large-spin-magnitude" 
systems mentioned above include systems in which one 
unit has more than two degrees of freedom, such as or- 
bitally degenerate systemsi^. In such systems, many 
coupling constants appear in general. However, exper- 
imental information about multipolar couplings is lim- 
ited. As a starting point to explore such systems, un- 
derstanding of one of the extreme limits must be use- 
ful. Therefore in this paper, we investigate an SU(n) 
symmetric caso 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 . 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 . 20 ! 21 . 22 in 
three dimensions, where n is the total number of internal 
degrees of freedom, by means of high-temperature series 
expansion (HTSE). 

We consider a simple cubic lattice. Let each site take 
one of the n colors denoted by \a) with a — 1, 2, • • • , n. 
Using the Hubbard operator X a(3 := |a)(/3|, the exchange 
operator is expressed as 



(i) When n = 2, this model is reduced to the ordinary 
Heisenberg model with s = 1/2 by relation 2Pij — 1 = 
4s; • Sj. 

(ii) The SU(3) case corresponds to s = 1, competing 
quadratic and biquadratic exchange interaction, 7 



1 + Py = s; • Sj + («, • sj) 5 



(3) 



(iii) The SU(4) Heisenberg model is related 
to spin 3/2 systems but more often discussed 
in the context of orbital- and spin-degenerate 
systema 10 . 11 ! 12 . 13 ! 14 . 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 . In particular, 
the model in three dimensions has been used as an effec- 
tive model of CeBg to explain magnetic-field dependence 
of the transition temperature of an antiferro-orbital 
orderingi22i2i The four local states can be represented by 
l+)H. Hl+L Hl-L where |±) and |±] repre- 
sent an orbital state and a spin state, respectively. The 
pseudo-spin operators are defined by t z \±) = ±||±), 
«±|T) = l±), s*\±] = ±i|±], a±|T] = |±], and the 
exchange operator for n = 4 is rewritten as 



PiJ = U ■ tj + Si • Sj + 4(*i • *j)(s; • Sj 



(4) 



The Pauli matrices r 
expression, i.e., 



2t and cr — 2s may simplify this 



4Py = Ti • Tj + CT; • <Tj + (Ti • Tj)(<Ti • <Xj) + 1. 



(5) 
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X 



a=l 0=1 



(1) 



Colors of sites i and j are exchanged when Py is applied. 
Then, an SU(n) symmetric Hamiltonian reads 



H:= J^Py, 

<iJ> 



(2) 



where the summation is taken over all the nearest neigh- 
bor pairs. We consider the antiferromagnetic case, J > 0. 
Let us show the relations with spin operators explicitly 
for some of the special cases below. 



Note that there is a different representation of SU(n) 
Heisenberg model studied in detail using Quantum 
Monte Carlo method (QMC)24. However, the QMC for 
the Hamiltonian J5J suffers from minus sign problems in 
more than one dimensionsii24 

If the number of competing order parameters is large 
and frustration exists, the transition temperature can 
greatly be reduced from the mean-field value even in 
three dimensions. That will be the case with the Hamil- 
tonian J5J). First, it is isotropic with respect to n 2 — 1 in- 
dependent interacting components. Furthermore, it con- 
tains frustration as most clearly seen in Eq. (J5J of the 
SU(4) case. Each of the 15 components, r°, a 13 , t 7 ct 5 , 
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attempts an antiparallel alignment. However, it cannot 
be attained simultaneously, e.g., simultaneous Neel states 
of t z and cr z produce a ferromagnetic alignment of the 
product r z a z . Such frustration is not so explicit in the 
SU(3) case in Eq. @ as in the SU(4) case, yet appear- 
ance of the minus-sign problem in the QMC indicates 
that similar frustration lies in it^a 

This frustration should become stronger as n increases. 
One can see in one dimension how the system finds a com- 
promise against this frustration. The SU(n) Heisenberg 
model in one dimension can be exactly solvedi^SiLSiS 
An important point is that the ground state has a quasi- 
n-site periodicity. The SU(4) case is numerically stud- 
ied in more detail; ground-state properties by the den- 
sity matrix renormalization group 10 , and thermodynamic 
properties by the QMGii and the HTSE±2J£. The real- 
space spin-spin correlation function as a function of spin- 
spin distance has a positive value every four sites. Its 
Fourier transform has a outstanding cusp at q = n/2 
and a small cusp at q = n. Regarding the temperature 
dependence, as temperature decreases correlation with 
two-site periodicity develops, and then at lower temper- 
ature another correlation occurs with n-site periodicity. 12 

Such peculiar correlation in one dimension could sug- 
gest an exotic ordering in higher dimensions. However, 
not much is known about the antiferromagnetic SU(n) 
Heisenberg model in three dimensions. In fact, its fer- 
romagnetic variant is studied in Ref. |2^ by the HTSE 
for the uniform susceptibility that needs less effort to 
be calculated than susceptibility of other wave- numbers. 
To our knowledge, our calculation in this paper is the 
first HTSE aiming at antiferromagnetic exchange of this 
model. In this study, we investigate the model system- 
atically by changing parameter n of SU(n), and report 
some unique features of spatial correlation at finite tem- 
perature. Namely, in Sec.[H]we show that the behavior of 
correlation functions is similar to that in one dimension, 
and in Sec. IIIII that the Neel temperature disappears as 
n increases from n = 2. 



II. TEMPERATURE DEPENDENCE OF 
CORRELATION FUNCTIONS 

What should be calculated here is (-Py), which is re- 
lated to a correlation function. For example, when n = 4, 
fifteen components in Eq. J^J contribute equally, and thus 

W) = <°f of) = Wofrtf) = (4Py - 1)/15. For 
general n, we can define correlation function between 
sites i and j as 



1 



1 



'A.. 



1 



(6) 



for i 7^ j,a ^ (3, which does not depend on a nor (3 
because of the SU(n) symmetry. Its Fourier transform is 
denoted by 5(q). 

The high-temperature expansion is performed by ex- 
panding the Boltzmann factor in (3. In practice, 




FIG. 1: Temperature dependence of the next-nearest- 
neighbor correlation function, 5i_j with i — j = (110). 



the series coefficients in the thermodynamic limit are ex- 
actly obtained by a linked-cluster expansion^ To obtain 
the series expansion of (Pij) up to 0[((3J) M ], we need to 
calculate Tr[(H L ) m ] and Ti[P id {U^) m ] for < m < M, 
where 7Yl is the Hamiltonian in a linked cluster. In calcu- 
lating the traces, we use a property of the permutation 
operator, and calculation of the traces is reduced to a 
combinatorial problem to count the number of circular 
permutations in a product of permutations>i2i22i 2 & Here, 
an important point of our analysis is that the series co- 
efficients are obtained as polynomials of n, for example, 



S(7r,7r,7r) 



n n l n 6 



(216- 22 n 2 



-{(3Jf 



(7) 



Namely, the order of the series for every n is the same. 
We have obtained the series for S*i_j up to 0[(f3J) 9 ] for 
arbitrary i— j, and consequently 5(q) up to 0[(/3J) 9 ] for 
arbitrary q. As a special interest, S(n,Tr,ir) is obtained 
up to 0[(pJ) 10 }. The series are extrapolated using the 
Pade approximation (PA). 

First of all, in order to see the temperature-dependent 
nature of the spatial correlation, we analyze the next- 
nearest- neighbor correlation function Shq as we have 
done in one dimension in Ref. Figure ^ shows the 
results. Both the axes are scaled so that the high- 
temperature limit of every n matches. Here, we have 
simultaneously plotted extrapolation from several differ- 
ent choices of the PA, and the difference between them 
approximately represents an error of the extrapolation. 

The lowest order of the series of Si_j has a Neel-order- 
type correlation in any n. That is, the series of Si x i i z 
starts withOpJ^I+liyl+H*!] with sign (-1)^1+^1+1^1. 
Therefore, Shq > at high temperature. However, as 
antiferromagnetic correlation of each interacting compo- 
nent becomes larger, each short-range order disturbs an- 
other because of the frustration. The change of the sign 
of Si io suggests that the correlation acquire a longer pe- 
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riod at low temperature. For n > 6, it is clear that Suo 
changes the sign at a low temperature that increases with 
n in this scale. For smaller n, the relevant temperature 
range is below the converged region, and it is difficult to 
conclude from this data. 

Such a change of correlation at low temperature also 
appears in the Fourier transform of the correlation func- 
tion. A naive extrapolation of the series for ^(q) shows 
bad convergence for several q. This is probably because 
information at position x is lacking when cos(x • q) ~ 0. 
In order to avoid it, we extrapolate the series of a 
complex function ^2 X>0 (X? ^ Xj" x )e~ lx -' cl w(x.), and af- 
ter that take the real part of the extrapolated function. 
Here, the summation is taken for x > 0, y > 0, z > 0, 
and w;(x) is two to the power of the number of nonzeros 
in x,y,z. This method also makes use of its imaginary 
part, and the convergence becomes better than the naive 
extrapolation. 

As we have seen in the next-nearest-neighbor correla- 
tion function, the analysis becomes easier as n increases. 
Therefore, we show 5(q) of the SU(16) case in Fig.[21for 
a diagonal direction in the q-space, S(q,q,q). Here, its 
high temperature limit, which does not depend on q, is 
subtracted. As temperature decreases, antiferromagnetic 
correlation develops and S(tt, 7t, tt) increases as the curve 
at T = 2J shows. However, with further decrease of 
temperature, S(7r, 7r, 7r) starts decreasing and 5(q) with 
other q increases. It is not very clear from Fig. if the 
maximum of S(q) starts moving from (tt, tt, tt). Hence, 
let us show another quantity. Note that if the second 
derivative of S(q) at q = (tt, tt, tt) changes sign, the po- 
sition of the maximum clearly moves. Its temperature 
dependence for the SU(16) model is plotted in Fig. [31 It 
shows a change of the sign around T ~ J . In fact, this 
temperature of changing sign seems to hardly depend on 
n of SU(n). The result above suggests that the Neel or- 
der disappear at least in the SU(n) model with large n. 
Namely, order with wave number q ^ (tt, tt, tt), or disor- 
der, should appear. Then, the next question is, at which 
n the Neel order disappears. 



III. NEEL TEMPERATURE 

It is known that the SU(2) Heisenberg model has the 
Neel order at low temperature. Therefore in this sec- 
tion, we gradually increases n of SU(n) from n = 2. 
As a preparation for that, we find a reliable way to es- 
timate a transition temperature first, using the SU(2) 
model. The Neel temperature Tn can be characterized 
by divergence of S(tt,tt,tt), namely, by a singularity of 
S(tt, tt, tt) as a function of (3 J. In order to analyze such 
a singularity, we use so-called D-log-Pade approximation 
(DLPA), i.e., the PA for the logarithmic derivative. In 
using the DLPA, transformation of the expansion vari- 
able may improve the convergence of extrapolation. We 



S(q)-S(q)| T 
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FIG. 2: The Fourier transform of the correlation function 
as a difference from its high-temperature limit, S(q, q, q) — 
S(q)|T-oo, for the SU(16) model. 
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FIG. 3: Temperature dependence of d 2 S(q, q, q)/dq 2 at q = tt 
for the SU(16) model by two different PAs. 



choose a transformation2L2£ 

(3J 



1 — a 2 x 2 ' 



(8) 



where a is an adjustable parameter that improves con- 
vergence. In fact, the singularity closest to the origin in 
the original series is near the imaginary axis, while Tn 
corresponds to a singularity on the real axis. With this 
transformation, singularities on the real axis approach to 
the origin and those near the imaginary axis go away. 
Since the DLPA can estimate a position of the nearest 
singularity most accurately, errors of the DLPA can be- 
come smaller by this transformation. 

In order to find an optimal a, we calculate Xn as a 
function of a for a couple of choices of the DLPA, and we 
adopt a at which difference among different DLPAs is the 
smallest. Figure 01 shows Tn as a function of a for three 
different choices of the DLPA. Here, [m/n] denotes the 
DLPA with a polynomial of order m over a polynomial 
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FIG. 4: The Neel Temperature T N of the SU(2) Heisenberg 
model obtained by three different DLPAs as a function of the 
parameter a. 
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FIG. 5: The Neel Temperature T N of the SU(n) Heisenberg 
model as a function of n with a = 1.04. Here, [5/4] and [4/5] 
are plotted together. 



of order n. Since [4/4] is from the series one-order lower 
than the others, we choose a at which [5/4] and [4/5] are 
the closest, namely, a = 1.04. Then, Tn/J obtained here 
is 1.898, which is close to those in the literature, 1.892 
by the QMC 29 and 1.888 by the HTSE 30 . In addition, 
we have obtained a critical exponent simultaneously. We 
assume that the critical exponent does not depend on the 
spin magnitude, and compare it with that of the classical 
antiferromagnetic Heisenberg model, in which S (tt , tt , n) 
is identical to the staggered susceptibility equivalent to 
the 'uniform susceptibility of the ferromagnetic model'. 
Hence we can compare the critical exponent 7. The es- 
timation from our calculation above is 7 = 1.399, which 
is close to 1.396 by a Monte Carlo method^, 1.406 by 
HTSBS, and 1.388 by a field theorjiS, of the classical 
0(3) model. Therefore, we trust this way of analyzing 
Tn, and use it also for the SU(n) model with n > 2. 
Also for n > 2, the variable transformation JSJ should 
work because the singularity closest to the origin in the 
original series is near the imaginary axis. 

Originally n is an integer because it is the number of 



internal degrees of freedom. However, after obtaining 
the series as an analytic function of n, we can regard 
n as a continuous variable that has a physical meaning 
only when it happens to take integer values. In fact, as n 
changes continuously, the properties of the series change 
also continuously. Therefore, we gradually increases n 
from n = 2 and see n-dependence of Tn. Figure [S] shows 
the results. Here, we fix a = 1.04, and we plot [5/4] and 
[4/5] together. The optimal a may depend on n. How- 
ever, since the difference between the two curves is very 
small, we regard it as an optimal a. As n increases, Tn 
decreases almost linearly, and at n ~ 2.45 the singularity 
corresponding to Tn runs away from the real axis to a 
complex value with a finite imaginary part. For larger 
n, we do not find any singularity on the antiferromag- 
netic side of the real axis. Therefore, we conclude that 
the value of n at which Neel order disappears lies in the 
range 2 < n < 3. 

In addition, we have also analyzed S(q) with q =/= 
(jr, 7T, 7t). However, we have not found any symptom of 
ordering with q ^ (tt, it, tt) in the temperature region 
that we can reach by the present order of the series. 



IV. SUMMARY 

In summary, we have performed high-temperature se- 
ries expansions for the SU(n) Heisenberg model in three 
dimensions with arbitrary n. First of all, we have cal- 
culated a nerrf-nearest-neighbor correlation function. At 
least at large n, it changes the sign at low temperature, 
which suggests that the correlation should not be like 
Neel order, but have a longer period. Analysis of the 
Fourier transform of the correlation function also sup- 
ports that the ground state of the large-n SU(n) Heisen- 
berg model does not have the Neel-order correlation. 
Then, we have turned to an approach from n — 2. Since 
the SU(2) Heisenberg model has the Neel order, it should 
disappear at a certain n. We have first found a reliable 
way to estimate a transition temperature by analyzing 
the divergence of the (7r,7r,7r) component of the corre- 
lation function. Next, we generalize n to a continuous 
variable and increase n gradually from n = 2. We have 
concluded that the Neel ordering disappears for n > 2. 
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